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i SUMMARY 

! 

j 

| The use of the NASTRAN program in the Shuttle wing stability analysis is 
.described. Details of the actual structure, the finite element idealization, 

^nd the NASTRAN results are given. A comparison of the NASTRAN results with 
( those obtained with another computer program and with hand generated results 
indicates good agreement. An alternate approach for solving eigenvalue problems 
is illustrated and shows a considerable savings in computer time. Some em- 
phasis is placed on the relationship of the NASTRAN analysis in the design pro- 
cess bringing out more clearly the contribution of the results and showing the 
importance of the mode plots. Finally, a deficiency in the NASTRAN plate ele- 
ments when used to model structures made up of intersecting plates is discussed. 
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INTRODUCTION 


The Space Shuttle wing structure (illustrated in Figure l) is a low aspect 
ratio double delta with relatively light spanwise compressive load intensity (up. 
to 875 kN/m (5000 lbs/in.) ultimate). The external skin surfaces are required to 
be buckle free during vehicle operation in order to insure the integrity of the j 
Thermal Protection System (TPS) tiles. These tiles , which are composed of silic 
fibers, protect the vehicle's aluminum primary structure on reentry such that ; 
the maximum temperature is limited to 176°C (350 F). The design consists of 
closely spaced lightweight stringers which meet the dual requirements of a low 
cover weight and a buckle free surface. This is fairly typical of structures . I 
that have been built around the concept of the optimum compression panel, modi- 
fied for shear and lateral loading. The stringer that was selected is a roll- : 
formed stretched aluminum and is illustrated in Figure 2. The unique shape of 
the stringer (eight internal corners and only two hinged free elements) per- 
mitted the use of gages as low as 0 . 65 mm ( • 022 in.) while maintaining a test 
crippling strength of 379 MPa (55 ksi). 

Large chordwise loads (normal to the stringer) were identified when the 
thermal loads, resulting from temperature gradients in the primary structure due 
to reentry heating were calculated from a comprehensive wing-fuselage finite ele 
ment analysis. These loads which are approximately 175 kN/m (1000 lbs/in.) be- 
came critical for design. Modifications to the structure were thus required to 
overcome the general instability problem of the stiffened cover panels subjected 
to large chordwise loads. These modifications included the addition of a light 
riblet midway between the main ribs plus stiffening of the stringer cross sec- 
tion with bulkheads to prevent their distortion. 

This paper illustrates how the NASTRAN program significantly contributed to 
the efforts of the stress analysts in identifying the mode shapes and the criti- 
cal parameters controlling the buckling strength of the cover panels. 
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cr critical 

p plate 

SB stringer bending 

ST stringer torsion 
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DISCUSSION 


Structural Model 


A typical cover panel, spanning two adjacent spars and two adjacent ribs, 
shown in Figure 3. One quarter of the panel (ABCD) which included seven and on 
half stiffeners was modeled. Conditions of simple support were imposed on the 
rib AB to simulate a half sine wave mode existing between ribs on the continuou 
cover. Simple support conditions were also imposed on spar AD. Symmetric or 
antisymmetric boundary conditions were assumed along the symmetry planes BC and 
CD. The finite element model consisted of approximately 570 nodes and 600 plat, 
elements (QUAD2). The modified hat section was modeled as a Y section as shown 
in Figure 3 with BAB elements added at the base of the Y to obtain the appro- , 
priate properties. Although this model represented the basic stiffness of the 
section, it did not account for the local distortion which proved to be signifi- 
cant. Local distortion of the real stringer can drop the effective torsional ] 
stiffness GJ to 12# of the St. Vennant value. Connecting the two legs of the 
stiffener to form a Y produces a closed cell that has a GJ of about 50# of the 
theoretical value of the modified hat. Note that a Y section when modeled as 
shown in Figure 5 is virtually distortionless due to the built in truss feature 
and hence the full torsional stiffness calculated by the St. Vennant equations 
for closed cells is appropriate. This closed cell value, however, is close to 
an effective value of 40# that would be obtained if the waist of the modified 
hat is included with the upper cell. This is consistent with the insertion of 
"blocks" to form stiffening bulkheads in the real stringer. 

Figure 4 shows a CALCOMP plot of the NASTRAN model where the cover skin ha 
been separated from the stiffeners for clarity. 

NASTRAN Analysis and Results for the Basic Structure 

For the NASTRAN analysis Rigid Format 5 was used, employing the inverse 
power method to obtain the eigenvalues and eigenvectors. 

The NASTRAN results, for the basic stiffened panel without modifications, 
are shown in Figures 5, 6, 7 and 8. Chordwise compressive loading (perpendic- 
ular to the stringers) was applied and two different sets of boundary condition 
were considered. For the modes depicted in Figures 5 and 6, symmetric conditio 
were enforced along boundaries BC and CD, while for the modes shown in Figure 7 
antisymmetric conditions were applied along BC and symmetric along CD. The fir 
buckling mode (for symmetric conditions) is shown in Figure 5 where the displac 
ments of the skin and the stringers are depicted. Also included is an end vie'v, 
at CD which clearly shows the number of half waves (m) and the location of node 
lines. Figure 6 shows the next four modes with symmetric conditions. It is nc 
that the first four modes are coupled; that is, they involve both bending and 
torsion of the stiffeners while the fifth mode is a pure bending mode since the 
stiffeners bend only. In a similar fashion the first five buckling modes for t 
antisymmetric conditions are illustrated in Figure 7. The NASTRAN results have 
been plotted in nondimensional form against m, the number of half waves in the 
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ph of the panel, in Figure 8. Note that while the results for the symmetric 
| antisymmetric boundary conditions are very close, they do not overlap. The 
son for this is that the stresses in the prebuckled state are slightly differ- 
i due to the difference in the boundary conditions. 

Comparison of NASTRAN Results 

The NASTRAN results can be compared with results obtained with another com- 
er program, VIPASA, which is described in Reference 1. This program treats 
smatic structures consisting of long plates connected along longitudinal lines, 
s the cover sheet and stringers were modeled with a series of 0.762-m- (30 in.) 
g, flat plates. The VIPASA model did not employ a Y-section representation 
the stringer, as did the NASTRAN model, but followed closely the actual geo- 
ric outline accounting for chem milling patterns as shown in Figure 2. The 
st five VIPASA mode shapes, along with the critical loads, are shown in Figure 
|for the case of symmetric boundary conditions. The mode shapes obtained with 
VIPASA program are identical to those obtained with NASTRAN j however, the 
jer of the modes and the values of the critical loads differ. 

A further comparison of the NASTRAN results involved the use of hand con- 
ations. Two different approaches were employed. In the first, the stiffened 
L er panel was treated as an orthotropic plate, and in the second, the energy 
hod was used to obtain critical loads. Figure 10 illustrates the application 
orthotropic plate theory to the stiffened panel. The governing equations 
[e Reference 2) together with the assumed mode shape, which takes on a half 
le wave along the length of the stiffeners and varies in the other direction, 

! shown. Note that the orthotropic plate equations in Figure 10 neglect the 
[■sional stiffness of the stiffeners and coupling between the membrane and 
Lding stresses due to the shift in the neutral axis. The coupling effect can 
taken into account by employing a more general orthotropic plate theory as in 
i“erence 3- However, this leads to a complicated set of equations which must be 
.ved using a computer. Furthermore, in the calculations employed herein, v x 
5 taken to be zero and the stringer contribution to (Gl) (the plate torsional 

;idity) was neglected. Results are shown in Figure 8. Note that while the 
jrall trend agrees with the NASTRAN results in general, lower buckling loads 
; obtained with the difference decreasing with m. 

! 

i The use of the energy method is shown in Figure 11. Here a mode shape was 
feumed, and the strain energy due to plate bending, stringer uending and torsion 
pL the potential of the applied loads was computed. Application of the princi- 
b of minimum potential energy then resulted in a value for the buckling load, 
b'e different modes were assumed which are depicted in Figure 12 together with 
b respective buckling loads. Note that in this approach the interaction be- 
sen the stringers and the plate is not fully accounted for. As shown in 
ference 2 this interaction is a function of the mode shape. An approximate 
y of accounting for this effect is suggested by Timoshenko in Reference 4 and 
volves computing the moment of inertia of the stringers about the bottom of 
e plate. The torsional constant for the stiffeners was taken as 30$ of the 
eoretical value of the modified hat section so that consistency with the finite 
ernent model is maintained. Results obtained with the energy method are also 
otted in Figure 8 and show excellent agreement with the NASTRAN results. 
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Use of the ALARM Subprogram 


The aforementioned NASTRAN analyses had between 2724 and 2853 degrees of 
freedom. Computer time is listed in Table 1. Also of considerable importance 
the elapsed time or wall time which ranged from 5 to 10 hours for these runs, i 
With such extensive computer residence periods, the computing system reliabilit 
becomes an important factor. In fact several computer system "crashes" were ex 
perienced during this effort which resulted in lost calendar time. 

For these reasons, and in anticipation of analyses with an even greater j 
number of degrees of freedom, an alternate approach was investigated for solvin 
large eigenvalue problems. The concept was to use a very fast eigenvalue solu- 
tion algorithm developed by Ojalvo and Newman (Reference 5) which was implement 
into a working program for a NASA contract (Reference 6). This subprogram, cal! 
ALARM (Automatic Large Reduction of Matrices) employs an automatic tridiagonal 
reduction technique which is identical to the FEER routine described in Referen 1 
7. The procedure was to generate the stiffness and incremental stiffness matri 
in NASTRAN and then to use ALARM to obtain the desired number of eigenvalues an 
eigenvectors. The modes were then passed back to NASTRAN to be plotted. Table 
shows a comparison of the eigenvalues obtained from NASTRAN and ALARM. The 
savings in computer time is evident from the comparison depicted in Table 5 , wn: 
shows that the NASTRAN-ALARM-NASTRAN procedure reduced the computer cost by two : 
thirds. It should be noted that the ALARM values given in Table 2 are for the 
17th iteration, and that after the first iteration the first five eigenvalues 
agree with those in Table 2 to four figures. Thus, the number of iterations 
could be cut to 2 or 5 which would result in a decrease of one-half in computer 
time for the ALARM step. 
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THE RELATIONSHIP OF THE NASTRAN ANALYSIS 


TO SHUTTLE WING DESIGN CONSIDERATIONS 

In the previous discussion comparisons have been made between the NASTRAN 
psults and values obtained by using other approaches. It is advantageous to 
ptempt to tie the various analytical and design activities together to obtain 
more comprehensive view of the problem. Figure 13 illustrates the relationship 
f several activities that have taken place relative to the Shuttle Wing insta- 
Llity problem. 

Initially all stability calculations were performed using hand methods 
Energy techniques using assumed mode shapes). Obviously these methods are 
united, the results being only as good as the assumed mode shape. Attempting 
p guess at the lowest mode for overall instability can be difficult if not 
jnpossible (needless to say it is dangerous as well). Thus the NASTRAN analysis 
jas undertaken to provide an overall check on the assumptions that were employed 
a the hand calculations. As such, the first NASTRAN analysis duplicated the de- 
ign without the intermediate riblet, while the second analysis incorporated the 
iblet. The main purpose of the riblet is to prevent the mode that consists of 
lternate bending of the stringers. Having prevented this type of motion, it was 
iscovered from the NASTRAN analysis that the stringer was still twisting as it 
assed over the intermediate riblet. (See Figure 14(a).) In order to prevent 
his, it was decided to control the clearance between the stud and the flange 
f the riblet (see Figure 15), such that when the stringer twisted it would 
ome into contact with the flange of the riblet. The riblet flange was modeled 
s a tie bar and the results of this modification are shown in Table 4 and 
igure 14(b) . 

I 

Cross-sectional NASTRAN plots indicated that the stringer section was still 
.istorting as it passed over the riblet. At this point we began to question the 
.dequacy of the NASTRAN representation of the torsional stiffness of the stringer 
Lnd in fact the effective GJ of the real stringer. To answer some of these ques- 
■ions a series of stiffener studies was initiated. An individual stringer and a 
102 m (4 in.) width of skin was modeled as an assemblage of bars, beams, torque 
;ubes and shear panels. This model more accurately preserved the geometric shape 
)f the real stringer. The model was analyzed for a sine wave distribution of 
torque applied at the skin stringer attachment and for a concentrated torque 
Applied at the end of the .38lm (15 in.) length. Results of this latter analysis 
J-ere correlated with test data. A cross-sectional plot of the model and the de- 
formed shape at different stringer stations is shown in Figure l6. Note the 
Revere bending distortion of the cross section at the loaded end and that the 
listortion decreases with the stringer station. If the stringer were long enough 
bhe distortion would completely disappear, an indication that the total torque 
rould then be taken by St. Vennant torsion. The sine wave torque loading (com- 
parable to the distribution of torsion that the stringer would have to resist to 
prevent buckling) was used to calculate an effective GJ value for the stringer. 
The calculated value was approximately 12 % of the St. Vennant torsional stiffness 
At this point several design changes were considered. These included adding 
"blocks" at various spacings which tied the two sides of the stringer together 
it the waist (see Figure 17(a)). This scheme tends to mobilize the upper portion 
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of the stringer. A second concept was to insert a vertical "spike" through the j 
stringer to tie the cover to the bottom and two sides of the stringer (see 
Figure 17(b)). This scheme is far more effective than the block scheme but alsol 
weighs more. A third possibility was also initiated, which was to investigate 
an alternate stringer concept that would not be subject to large distortion undel 
torsional loading. 

The Shuttle Wing design 101 used a combination of blocks and spikes. A 
spike was placed adjacent to each of the riblet attachment points. This con- 
figuration was investigated by modifying the previous NASTRAN model. The result] 
of the analysis are shown in Table 4. The final NASTRAN analysis consisted of 
adding a beam element to the skin parallel to the intermediate riblet. This las^ 
modification which prevented bending of the skin along the riblet did little to 
increase the buckling allowable. 

Some concern was developed at this time over the adequacy of elements in 
NASTRAN to predict the buckling behavior of structures composed of QUAD2 plate ! 
elements that intersect at an angle. In the Q.UAD2 elements, the inplane motion 
is assumed to be linear between node points while the out of plane motion is 
essentially a cubic. An incompatibility problem thus inherently exists such 1 
that elements connected at angles to each other do not continuously reinforce 
one another, but instead are "stitched" together at the nodes. In addition to 
this, the elements lack torsional stiffeners normal to their plane. Thus, as a 
separate study, the QUAD2 elements were used to obtain the buckling modes of an 
equal leg angle column. It was discovered that while the torsional buckling 
mode of the angle could be obtained using a refined grid, convergence was from 
the low side. Furthermore, some peculiar modes were observed and the Euler 
buckling mode could not be obtained. In order to clarify this problem, an in- 
vestigation to determine the applicability of higher order plate elements (as 
described in Reference 8) has been initiated. These concerns caused us to turn 
to an alternate analytical approach, the VIPASA program. VTPASA can model con- 
tinuous structure in detail. It cannot, however, account for discrete stiffen- 
ing or items that do not fit into a Fourier series expansion. The correlation 
of the VIPASA results has been discussed previously. 

As discussed above, work was initiated to develop an alternate stringer con 
cept with the prime goal of minimizing the local distortion. This involved the 
analysis of a typical unit width of stringer with a torsion load applied at the 
skin line which is balanced with a shear flow distribution. The stringer 

behaves as a frame, and hence has an internal force and moment distribution. 
Dividing the moment by the force gives the location of the load line. The geo- 
metric shape of the stringer was then modified to straddle the load line, keepip 
in mind other considerations such as local crippling requirements, fabrication 
restrictions and clearances dictated by the present design. The process of 
analysis and geometric modification was iterated until a stringer of the shape 
shown in Figure 17(c) was obtained. The effective GJ of this stringer is approx 
mately 5 times the effective GJ of the unreinforced modified hat. Work is not 
yet complete on this concept but strong consideration is being given to it for 
inclusion in future Shuttle Wing designs. 

A useful plot that illustrates the effects of the various design parameters 
is shown in Figure 18 where the critical stress is plotted as a function of the 
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lsupported stringer length. The plot was obtained using the energy equations 
.ven in Figure 11. The plot illustrates the improvement in allowable buckling 
-ress in using the frame stiffened hat over the modified hat. Buckling stresses 
re indicated for a range of modified hat effective GJ values; the particular 
ilue within the range depends on the combination of block and spike spacing and 
pringer length. Note that the effective GJ for a stringer varies with length 
id is equal to the St. Vennant value at L = » and virtually zero for L = O.It also 
lould be noted that from a practical point of view one can question the real 
: f ectiveness of the blocks and spikes as these quantities are related to various 
^nufacturing procedures. The four lowest NASTRAN values from Table 4 have been 
iperimposed on the plot at spacings of 0.762m (30 in.) and 0.38lm (l5in.). NASTRAN 
ialysis I did not include the riblet; analysis II did but did not prevent twist- 
ig at the riblet support (hence it is lower than the value predicted by the hand 
ilculations which tacitly assume zero rotation at the supports). NASTRAN analy- 
Ls III prevented rotation at the support but not distortion. Analysis IV pre- 
bnted rotation and distortion by inserting a spike at the riblet. 

[ The plot also indicates the improvement in buckling stress due to inserting 
ae riblet (spacing changes from 0.762m (30 in.) to 0.38lm (15 in.)) which prevents 
|ie first bending mode. 

It is important to note that simplified tests of the stiffened panel, under 
pans verse compressive loading, correlated well with some but not all of the 
lalytical results. In general, the large strength increases in the bending 
|id coupled modes due to the addition of the riblet were attained. Also, the 
nprovement in the effective torsional stiffness of the stringer with the block 
,id spike modifications was corroborated. For the torsion modes, however, an 
^fective length reduction from 0.762m (30 in.) to 0.38lm (15 in.) proved diffi- 
ilt to attain. Subsequent analysis, using hand computations , showed that these 
3wer than expected strength levels could be traced to the low torsional stiff- 
ess provided at mid bay. In effect the torsional restraint at the riblet was 
eo flexible and allowed twisting as the plate buckled. This agreed with test 
bservations. There results demonstrate the need to provide and model accurately 
le local stiffness which affect the behavior of the actual structure. 
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CONCLUDING REMARKS 


The use of the NASTRAN program in the stability analysis of the Shuttle . 
wing panels has been described. Moreover, the role of the NASTRAN analysis in 
the design process has been discussed to illustrate how it impacted the design. 
The NASTRAN results helped to identify the critical modes and show the impor- 
tance of the torsional properties of the stiffeners. In this effort, the NASTRAJj 
plots of the mode shapes were indispensible in that they indicated where the 
fixes should be made and clearly demonstrated their effectiveness. The results 
indicated that for future analyses (which will include spanwise loading in 
addition to the chordwise loading to determine the interaction effects) a better 
representation of the stringers is required. ; 

It was demonstrated that a tridiagonal reduction scheme can greatly reduce j 
the computer cost in solving large eigenvalue problems. | 

Finally, it was discovered that the plate bending elements in NASTRAN have : 
incompatibilities which can lead to difficulties when the elements are used to ' 
model structures made up of intersecting plates. Further work in this area is 
clearly required to develop a plate bending element which would overcome the 
present deficiencies. 
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igure 2. - Shuttle wing cover pane 



0.762m 








Mode 2- x = 1.666, N cr = 100. 5 kN/m (574 lb/in) m=9 



Mode 3 - X =1.154, N cr =109.2 kN/m (623 lb/ in} m=l; 




Mode 5 - X = 1.0, N cr = 129.O kN/m (737 lb/in} m=l! 


Figure 6.- NASTRAN results - modes 2-5 - symmetric conditions 



Mode 1 - X =1-5, N cr =96-5 kN/m (551 lb/in) m=10 


| 
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Mode 5 - X = 0.9375, N cr = 13^-9 kN/m (770 lb/in) m=l6 
Figure 7*- KASTRAN results - modes 1-5 - antisymmetric conditions. 
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Figure 10.- Orthotropic plate equations 
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Figure 11.- Energy approach. 
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Figure 12.- Results of energy approach. 
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Figure 13.- Relationship of design analysis activities 











NASTRAN ANALYSIS - CASE II 1 st MODE 
N cr = 146.2 kN/m (835 lb/in.) 
(A) 



NASTRAN ANALYSIS - CASE III 1 st MODE 
N cr m 215.2 kN/m (1229 lb/in) 

(B) 

Figure 14.- NASTRAN results for modified structure. 
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TORSION LOADING 



Figure 16.- Distortion of stringer. 
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Figure 18.- Effect of design parameters. 



